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We study thermodynamic properties of spatially separated electron-hole plasma in double-layered 
systems using Green function formalism. The screening of the Coulomb interaction is considered 
in the framework of Thomas-Fermi approximation, and a qualitatively new mechanism of screening 
by indirect excitons is taken into account. The exciton density is shown to decrease sharply with 
increasing electron-hole separation up to one exciton Born radius. The strong mutual enhancement 
of screening and charge-separation effects is found. 
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I. INTRODUCTION 

Spatially separated quasi-two-dimensional (2D) electron-hole plasma can be found in different modern nano- 
heterostructure-based devices. There is a continuous attention to closely situated (so called double) quantum wells. 
This attention is mostly due to the possibility of the Bose-Einstein condensation of spatially indirect excitons in 
such a system [0. Another physical object with spatially separated electron-hole plasma, which currently attracts 
increasing interest, is a GaN/AlGaN quantum well, grown on a sapphire substrate In the latter case electrons 
and holes are separated by huge (of the order of one MV/cm) internal electric field. Because of the enhanced binding 
energy in wide-gap semiconductors, excitons in nitride-based nanostructures are a subject of extensive experimental 
and theoretical research. 

The aim of this work is to provide a consistent many-body theory, describing statistical mechanics of the quasi-2D 
exciton/free carrier plasma with spatial charge separation. 

II. SCREENING IN THE SPATIALLY-SEPARATED ELECTRON-HOLE SYSTEM 

First of all, the potentials of electron-electron and electron-hole interaction should be established. The Thomas- 
Fermi (or, as it is often called in the 2D case, Stern-Howard) approximation is one of the most frequently used 
methods for dealing with static screening of Coulomb interaction in semiconductor systems [|J . In the frame of this 
approximation the charge, induced by external perturbation, originates from a redistribution of free carriers and is 
proportional to a local electric field. The resulting electric field potential is given by the Poisson equation in the form: 
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where <j) a {p) — J 4>{P: z )9a( z ) dz, g a (z) is a free-carrier-distribution function in the direction normal to the quantum 
well plane, k is a background dielectric constant, Q ex t is an external charge density, and the screening parameter in 
the 2D case is 
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Here n a is a free particle density, fi a is a quasi-Fermi level, m a is an effective mass, fi e h — m e mh/(jn e + rrih) is the 
reduced effective mass, as = nh 2 /([j,*e 2 ) is the 3D exciton Bohr radius, and (3 — l/(fcsT). 

The contribution of excitons to screening in the conventional single quantum well can be treated within the Thomas- 
Fermi approximation as a reduction of free carrier densities entering Eq. (|J) . In the case of spatially separated plasma, 
screening by excitons can not be neglected even in the Thomas-Fermi approximation. The indirect exciton has the 
finite dipole moment, and the external electric field causes redistribution of excitons, which plays its own role in the 



* Corresponding author; E-mail address: V.Nikolaev@exeter.ac.uk 



1 



resulting field reduction. Using the same approach as for a free-carrier plasma, the indirect-exciton part of the induced 
charge density can be written as 



Qtnd = ^ e s XC (9h(z) - g e (z)) [eMp) - e$ e {p)] ■ (3) 
Here we introduced the excitonic screening parameter 
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where n exc is the exciton density and /i = fi e + jih- 

In the present paper we consider a simple model of spatially separated electron-hole plasma, where electrons and 
holes are purely two-dimensional and reside in different planes with distance d between them. The generalization 
of this model to incorporate finite-width quantum wells will be presented somewhere else. The two-plane model is 
particularly valuable for studying the interplay between screening and spatial separation of different types of carriers. 
The problem of finding electron-electron repulsion and electron-hole attraction potentials has cylindrical symmetry, 
so we use the following representation: 
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where Jo is the Bessel function. Let us assume that the electronic charge resides in the plane with coordinate z = d 
and the holes are in the plane z = 0. In this model the distribution functions g a {z) are given by delta-functions. The 
Poisson equation is reduced to a one-dimensional problem: 
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Solving Eq.(|), we obtain repulsion and attraction potentials in ^-representation: 
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where both terms 3 = [1 - exp(-qd)] [2k + (q e s + q' s l )(l + exp(-qd))}, and 6 = (q + q { s e) )(q + q [ s h) ) - q { s e) q { s h) exp(-2gd) 
in denominators of Eqs.(0,||) do not depend on the excitonic screening parameter q e s 



~exc 



III. CORRELATED DENSITY 



In order to provide a consistent treatment of many-body effects in a 2D electron-hole plasma, we base our theory on 
a quasi-particle picture and Green functions technique. Electrons and holes are redistributed among the quasi-particlc 
states with a sharp spectral peak and dispersion e a = H 2 k 2 /2m a (where the energy is measured from the corresponding 
size-quantized level renormalised by the many-body shift) and some other "correlated" states originating from the 
interaction between quasi-particles. 

The total density of electrons can be divided into two parts, the first part representing "free" particles with a 
quasi-Fermi level fj, e , and the second part incorporating interaction between quasiparticles. This second part, which 
is called the correlated density, is given ^| by 

n{r = 2 £/ ^r e(fc ,.) [/e(M - fe (e e{m ^7M^- ^ 
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Here T e (k,uj) — ImS e (fc, w) is an imaginary part of the single-particle self-energy, f a {e) — [exp(/3(e — fx a )) + 1] 1 is 
the Fermi distribution function, and the derivative of the principal value is given 
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From the point of view of diagrammatic technique, exciton is a succession of electron-hole scatterings, which is repre- 
sented by the series of diagrams serving as a basis for the ladder approximation Q] . We use the ladder approximation 
to our many-body problem, keeping only ladder-type diagrams in the expression for the electron self energy: 
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where the summation is taken over different types of particles, b — e,h, 2D wave-vector k' and Fermi-type Matsubara 
frequencies Q b . The quantity T represents a sum of ladder diagrams for electron-electron or electron-hole scattering 
and is given by the frequency-dependent T-matrix equation: 
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where the potentials V^ b are given by Eqs.(0,||), 
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and 



N ab (k, q) 
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N ab (k, q) = l- f a (e a (k)) - f b (e b (q- k)). 

Following Zimmermann |^], who derived a similar formula for the 3D case, we obtain the expression for the correlated 
density in the 2D electron-hole system: 
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Here e mn are the values of the exciton bound state energies, where the first index m denotes the angular momentum 
and the second index n numerates different states for a given m. The bound state energies correspond to the poles of 
the T-matrix. The quantities S^(k) give the phases of the coefficients in the Fourier series expansion of the T-matrix. 
The values of these quantities coincide with the scattering phase shifts when the k-space filling is neglected. Thus, the 
second part of Eq.(|l2|) describes the scattering states contribution to the correlated density. In the low density/high 
temperature limit this expression differs from that of Portnoi and Galbright only by the factor 2 sin 2 in the 
integrand. The "e-e" part of the correlated density has a negative sign and is responsible for the reduction of the 
plasma density due to repulsion between the carriers of the same type. The electron-hole part of the density is positive 
and can be attributed to the presence of the excitons. As it was shown in Ref. Q for dilute plasma, this part of the 
density does not change abruptly when one of the bound states disappears with increasing screening. We treat the 
electron- hole part of the correlated density as the density of excitons: 
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Substituting n exc from Eq.(0) into Eq.ffl) we get the following expression for the excitonic screening parameter: 
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/s(e) = [exp(^(e-/Ze -///»))-!] ■ 

To calculate the exciton density we have to run the following self-consistent procedure. Firstly, for given electron- 
hole density we calculate the screening parameters from Eq.(||) and obtain corresponding bound state energies and 
phase shifts. Then, using these results we calculate the new values of excitonic screening parameter and free carrier 
density and recalculate bound state energies and phase shifts, repeating this procedure until the exciton density stops 
changing. 



IV. RESULTS AND DISCUSSION 



Throughout the calculations we neglect the phase shift dependence on the carrier density and obtain the values of 
S a b using the 2D modification |tJ of the variable phase approach. To find the bound state energies we use the method, 
proposed by Campi and co-workers || , which is also based on the ladder approximation. We take into account only 
the ground exciton state. The material parameters are chosen to match these of the GaAs/AlGaAs quantum wells, 
and the difference between dielectric constants of the QW and the barrier is neglected. 

In Figure 1 we present the dependence of the binding energy of the exciton ground state on the inter-plane distance. 
One can see (Fig. lb) that the decrease of the exciton binding energy with increasing the distance is substantially 
enhanced by the presence of free carriers and indirect excitons, which screen the Coulomb attraction. Figure 2 shows 
the exciton density as a function of the total plasma density for different inter-plane distances. One can see that the 
exciton density has a maximum and then decreases sharply with increasing the total number of carriers, which reflects 
the action of both screening and k-space filing. 
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FIG. 1. The ground state binding energy of the indirect exciton as a function of inter-plane separation for different 
free-carrier screening parameter: ql = (solid line), q% = 0.25 a^ 1 (long-dashed line), q% — 0.5 a^ 1 (dashed line), q% — a^ 1 
(dotted line), (a) - measured in the 3D exciton Rydberg units, (b) - normalized by its maximum value. 
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FIG. 2. Exciton density as a function of plasma density for different inter-plane distances: d = (solid line), d — 0.25 as 
(long-dashed line), d — 0.5 as (dashed line), d — as (dotted line). 
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